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INTRODUCTION 



Spatial distributions of objects are ubiquous in na- 
ture. The pattern of distribution of a population can 
be characterized by statistical relationships between each 
component and its nearby companions. Methods based 
on neighbor relationships and two-point correlations are 
some of the simplest and most popular techniques for 
statistical pattern recognition. The analysis of spatial 
patterns is especially useful in modeling and describing 
a wide variety of natural PHI] and social phenomena [J- 

Methods for spatial structure studies of many-body 
systems have been developed and extensively applied in 
physics of condensed matter 0, HH • Similar techniques 
are employed to understand the associations of galaxies 
and the large scale structures of the Universe [16|, Ht) ]. 
The main spatial structure function in the study of liq- 
uids is the pair distribution function (PDF). For statisti- 
cally homogeneous and isotropic systems, which are the 
focus of this paper, the PDF is a radial function g(r), 
traslationally invariant, which characterizes the proba- 
bility density of finding a particle at some given distance 
r from a reference particle. Assuming pairwise particle 
interactions, g(r) can be used to obtain equations in- 
volving macroscopic thermodynamic variables [l5j . An 
additional description can be obtained from /cth-nearest- 
ncighbor PDFs gk(r) [l8r- 22] . The first nearest neighbor 
distribution, g\ (r), is particularly useful to understand 
local fields around a particle in a fluid [23| - r26j . or a star 
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in a stellar cluster or galaxy 27 1, and it has applications 
in other problems in physics [281 ] . However, a further de- 
tailed microscopic information is sometimes required to 
statistically characterize physical properties of fluids and 
materials. For instance, mutual-ncarcst-ncighbor (MNN) 
pair descriptors are used in the study of solute relax- 
ation processes in fluids [2{|, phonon spectra analysis of 
disordered systems (30j . and for explaining the behavior 
of glass- forming liquids |3l| . In general, the probability 
that particles in a fluid form MNN pairs can not be ex- 
actly deduced from the standard g(r) or individual PDFs 
gk{r). Therefore, a superior spatial structure descriptor 
is highly desirable. 

We propose here the introduction of a more complete 
description of the microscopic neighborhood structure 
around a reference particle. We introduce a conditional 
pair distribution function (CPDF) which gives the prob- 
ability density of finding a particle to a certain distance 
from other in a fluid, with each of these two having a 
nearest neighbor to a fixed but otherwise arbitrary dis- 
tance. This function describes special four-body config- 
urations, but also contains contributions due to the mu- 
tual nearest neighbor (two-body) and shared neighbor 
(three-body) configurations. The new PDF is a spatial 
descriptor at the two-point distribution level but con- 
tains nontrivial higher-order structural information. In 
particular, it is a superior descriptor sensitive to topolog- 
ical connectedness information of the point pattern and 
could resolve a variety of degenerate point configurations 
associated with the standard PDF g(r) [32| . 

The primary purpose of introducing the CPDF is to 
establish a correct expression for evaluating the interpar- 
ticle potential energy in the so-called thermostatistical 
space partition (TSP) formalism for fluids. The TSP ap- 
proach, which was introduced in Ref. HH for the case of 
pure substances and extended to mixtures in Ref. [34| . 
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was motivated by difficulties with the calculation and 
qualitative prediction of optical properties in white dwarf 
stars, currently based on the occupation probability for- 
malism of Hummer & Mihalas 3^ and whose consistency 
has been discussed in Ref. [36 1. The TSP method com- 
bines the free-energy minimization technique (37j with 
space partitions that asign to each particle a volume v, 
which is determined by the distance of the particle to 
its closest neighbor, and roughly measures the free space 
surrounding each particle. The CPDF introduced in the 
present paper, denoted g vv >(r), is a PDF where the vol- 
umes v and v' of both pair particles are specified. It will 
be shown here that this CPDF is needed for an appropri- 
ate representation of the potential energy of fluids within 
the TSP formalism. 

Explicit evaluation of the CPDF g vv < (r) is considerably 
more complicated than the computation of the conven- 
tional PDF g(r). In the pairwise additive assumption for 
the potential energy of a fluid, the CPDF can be eval- 
uated using a hierarchy of integral equations for s-body 
correlation functions, for instance the Yvon-Born-Green 
hierachy fl5j , truncated by introducing some appropriate 
closure approximation. In the general case, with g vv > (r) 
representing four-particle configurations, the hierarchy 
expansion can be truncated in factors with correlation 
functions higher than s = 5 (in analogy with the mean- 
field approximation used in Ref. [29| to assess the MNN 
pair distribution of Lennard- Jones fluids) . A method of 
calculating g vv i (r) by means of this scheme is currently 
being prepare in detail. In the present contribution we 
will restrict ourselves to the evaluation of g vv > (r) for ran- 
dom, non-interacting systems or so-called Poisson ensem- 
bles, which can be directly derived from simple and well- 
known probabilistic properties (e.g., @,[3al39||)- The aim 
is two-fold. First, the study of the CPDF for random 
points may prove to be a useful guide for investigating 
this function in other systems, apart from its own impor- 
tance at a fundamental level. Second, Poisson statistics 
is an especially powerful tool for studies of a wide variety 
of systems, and its CPDF may provide a new method of 
analysis. In fact, Poisson statistics emerges as a useful 
alternative in view of the practical limitations of a rigor- 
ous treatment for the complex natural systems. Random 
samples are employed very widely by comparison or con- 
trast in the analysis of observed event data doj, or as 
a convenient first approximation to describe the spatial 
distribution of a variety of natural systems, from biologi- 
cal species (J-Q to stars on the sky (4l| - |43j . The present 
work is part of an effort dedicated to improving the TSP 
formalism and to providing a well defined and very de- 
manding test model of the CPDF for further research. 

The paper is structured as follows. In Section [TT1 the 
conditional pair distribution function is formulated and 
general results are presented. Its implication within the 
TSP formalim is described in Section Hill Section llVl con- 
tains a derivation of the CPDF for Poisson patterns in 
one, two, and three spatial dimensions. Some derived 
results arc examined there. Concluding remarks are con- 



tained in Section iVll 

II. DEFINITIONS AND GENERAL RELATIONS 

Consider a system of N identical particles distributed 
in a D-dimcnsional Euclidean space with volume V. For 
conciseness, the radial distance r from a reference particle 



will be expressed in terms of the volume u> 
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of a 



.D-dimcnsional sphere of radius r centered in the particle, 
with 



r D/2 



r(i + #)' 



(2.1) 



where T is the gamma function. The conventional PDF, 
henceforth expressed as g(u>), describes the probability 
density of finding a particle on the surface of a sphere of 
volume uj centered in the reference particle. More exactly, 
g(u>) is defined such that 

ng{ui)duj = mean number of particles between the 

surfaces of uj and uj + dui centered at the 
reference particle, (2-2) 

where n = N/V is the mean density in the fluid. 

It is assumed here that for any reference particle there 
is one and only one nearest neighbor. This means that 
if two or more particles are equidistant from the refer- 
ence, one of them will be arbitrarily designated as the 
first neighbor. The distance a between a particle and its 
nearest neighbor will be expressed by the volume v of a 
D-dimcnsional sphere of radius a. The volume v will be 
referred to as the available volume (AV) of the particle 
because it roughly measures the free-particle space sur- 
rounding it. The size distribution of AV in the fluid is 
given by the distribution n v , where n v dv is the number 
density of particles with available volume between v and 
v + dv. The quantity (n v /n)dv is clearly the probability 
that the nearest neighbor lies in the volume dv of a spher- 
ical shell centered in the reference particle. Therefore, 
n v /n is equivalent to the so-called first neighbor distri- 
bution (Udfl. In the thermodynamic limit (TV, V — > oo, 
keeping N/V constant), it follow that 



n v dv. 



(2.3) 



The conditional pair distribution functions g v {w) and 
g vv i{u>) are defined by the following quantities 

ng v (uj)duj = mean number of particles between the 

surfaces of uj and u) + duj centered at the 
reference particle with AV v, (2-4) 

n v i g vv i [ui)dv' dui = mean number of particles with AV 

between v' and v' + dv' lying between 
the spherical volumes uj and uj + dui 
centered at a particle with AV v. 




FIG. 1: (Color online) Left: Example of a point pattern with 
a characterization of the evaluation of the pair distribution 
functions g(to), g v (ui) and g vv i(uj). The point P' lies in the 
surface of the spherical volume ui centered at the reference 
point P. Nearest neighbors of P and P' lie in the surfaces of 
spherical volumes v and v' centered in P and P , respectively. 
Right: Representation of the four classes of (v, v') pairs (see 
text). Here, v+ = max(v,i/) and V- = min(v,v'). 



(2.5) 

According to the meaning of n v , it is clear that g(ui) is 
just a weighted average of g v (uj) over the whole range of 
V, that is 

/>oc 

ng(uj) = I n v g v (cu)dv. (2-6) 
Jo 

Similarly, the integration of (|2.5[) over v' covering the 
whole range of available volume gives (|2.4p . i.e., 

/•OO 

ng v (uj) = / n V 'g VV '(uj)dv'. (2.7) 
Jo 

From (|23| and (|2~7|) it follows that 

f>00 f 1 oo 

n 2 g(u)) = / n v n v rg vv r(ui)dvdv' . (2-8) 

Jo Jo 

From the definitions given above, we can see that the 
function g v (uj) is the probability density that a particle 
lies at the surface of a spherical volume to centered in a 
randomly chosen particle with AV v. Similarly, g vv >(uj) is 
the probability density of finding a particle with AV v' in 
the surface of a spherical volume uj centered in a particle 
with AV v. The conditional character of the probability 
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TABLE I: Classes of pairs (v,v') and some properties. 



Configuration 


n-body 


LI- 


range 


(i) mutual nearest neighbor 


2 bodies 


D 


/ 

= V = V 


(ii) simple nearest neighbor 


3 " 


D 


= v+ = max(w, v') 


(iii) shared nearest neighbor 


3 " 


<'H 


- < UJ < UJ vv i 


(iv) other 


4 " 


'<M 


- < UJ 



densities g v (u) and g vv /(uj) appear due to the specifica- 
tion of the AV v and v' for one or both particles of each 
pair. Due to the symmetry of the problem, the function 
g vv '(u>) is symmetric in the variables v and v' . Fig. [T] 
(left) displays a schematic point pattern with characteri- 
zations of the standard PDF g(u>) and the CPDFs g v (ui) 
and g vv > (w). The function g vv i (w) and the distribution n v 
are the central quantities of the present issue since knowl- 
edge of them allows one to calculate g v (u->) and g(u>) with 
(|27T]) and (im 

The function g vv / (uj) describes in general four-body 
configurations, i.e., pairs (v,v') composed of two parti- 
cles with a third particle in the surface of the volume v 
centered in the first one and a fourth particle in the sur- 
face of v' centered in the second one. However, there are 
also other configurations of pairs (v,v') formed by only 
two or three particles. As is illustrated in Fig. [1] (right) 
and detailed in Table[H there are four independent classes 
of pairs (v, v'): 

• Class (i) denotes mutual nearest neighbors, i.e., 
each particle of the pair (v,v') is the nearest neigh- 
bor of the other. 

• Class (ii) corresponds to the configuration where 
one and only one particle of the pair (v,v') is the 
nearest neighbor of the other. 

• Class (iii) refers to the shared nearest neighbor con- 
figuration, i.e., both particles of the pair (v,v') have 
the same nearest neighbor. 

• Class (iv) denotes pairs (v,v') which involve a con- 
figuration with four different particles. 

Because classes (i-iv) arc mutually exclusive and exhaus- 
tive, we have 

g m ,{u) = g»( U ) + + 9 ^(uj) + gW(u,), (2.9) 

where g^J,(uj) is the contribution to g vv i(uj) due to pairs 
(v,v') of class (v). 

Geometrical considerations which involves the avail- 
able volume v of a particle or, equally, the distance to its 
nearest neighbor, can be employed to derive some gen- 
eral properties of the contributions to g vv i (ui) on the right 
hand side of (|2.9[) . As can be seen from Fig. [TJ pairs class 
(i) only occurs for u> = v = v' and, therefore, the math- 
ematical representation of the term g^, (uj) must have 
singularities. Indeed, taking into account (|2.5p . it is eas- 
ily demonstrated that this term is proportional to the 
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product of two Dirac delta functions, 

9vl'( u ) « K v ~ v')8(u - v). 



(2.10) 



A similar behavior is deduced for the second term on the 
right-hand side of (|2.9[) . Pairs (v,v') of class (ii) occur 
only for to = v+, with 



and therefore 



v + = max(t>, v'), 



9vv'(u) oc(5(w -v+). 



(2.11) 



(2.12) 



Since two particles cannot share the nearest neighbor 
when their volumes v and v' do not overlap, the third 
term on the right of (|2.9jl vanishes for uj > w vv i , where 



(2.13) 



The third term vanishes also at uj < v+ because no parti- 
cle lies at lower distance than the nearest neighbor. Thus 



9vv' ] ( w ) = °> u i ( v + > w ™ ' ) ■ 



(2.14) 



Similarly, the fourth term vanishes for uj < v + , 

g^)(cj)=0, uj<v + . (2.15) 

Therefore, g vv > (uj) vanishes at uj < v + and has singulari- 
ties for pair configurations with uj = v = v' and oj = v+. 

At the fluid context and in analogy to the standard 
PDF, one may expect that, at sufficiently large interpar- 
ticle distances (large uj), pairs (v,v') are uncorrelated and 
the CPDF must then reduce to unity 



lim g VV '(uj) = 1. 



(2.16) 



III. THE ROLE PLAYED BY g vv ,(u) IN THE TSP 
FORMALISM 



Let us define u vv > (uj) as the interparticle potential be- 
tween two particles forming a pair (v, v'), with a particle 
centered in the volume uj and another one located on the 
surface of uj. The potential u vv >(oj), where the AV v and 
v' for both pair particle are specified, is a generalization 
of the usual pairwisc interparticle potential u(uj). u vv > (uj) 
formally takes into account perturbative effects on the 
two-body interaction potential introduced by the nearest 
neighbor of each particle in the pair. Many-body contri- 
butions to effective pair potentials are known for liquid 
hydrogen |44| - |46l| an d have been theoretically considered 
for fluid helium 47]. However, potential energy surfaces 
of four-particle systems are needed for a full knowledge 
of u vv i (uj) in a fluid. For the moment, ab initio potential 
evaluations are limited to three-body systems at few ge- 
ometrical configurations and for some atomic, molecular 
and colloidal systems [48| - |50} . 



Multiplication of (|2.5I) with u vv >(oj) and integration 
over all uj, v and v' , yields the total potential energy 
of the system 



2V 



V r V 



JO 



N V N V , 



(w)flw (oj)doj 



dv'dv, 



(3-1) 

where N v = n v V, N v > = n v >V, and the factor 1/2 avoids 
the double count of pairs (v,v'). If one writes 



u v (uj) 



J n v >u vv >(uj)g vv '(uj)dv' 



J n v >g vvl (uj)dv' 
then (|3.1[) can alternatively be written as 



(3.2) 



N 



u v (uj)g v (uj)duj 



dv. 



2V J 

Similarly, if we define an averaged pair potential as 

/ n v u v (uj)g v (uj)dv 



u(uj) 



J n v g v (uj)dv 



(3.3) 



(3.4) 



then the total potential expression reduces to that typ- 
ically used in physics of liquids, which is based on the 
conventional PDF g(uj), i.e., 

N 2 f 

$ = — / u(uj)g(uj)duj. (3.5) 



2V 

For systems with two-particle potential without many- 
body effects, u vv '(uj) = u v (uj) = u(uj). Otherwise, the 
CPDF g vv i (uj) can be used to calculate effective pair po- 
tentials [Eqs. (|3.2[) and (|3.4|) ] directly from ab initio 
many-body potentials u vv >(uj) (when these become avail- 
able), and thus the conventional evaluation of the total 
potential energy in liquid theory is preserved [Eq. (|3.5[) ]. 

The thermostatistical space partition (TSP) formal- 
ism is a mathematical model for moderately dilute gases 
which combines the free-energy minimization method 
[37j with space partitions. It was developed for one- 
component [33| and multi-component [3~i| fluids. Appli- 
cations to simple gas models (hard spheres, van der Waals 
fluids, attractive particles, and partially ionized hydrogen 
gas) have shown that TSP provides a useful tool to evalu- 
ate self-consistcntly thermodynamic and structural prop- 
erties of gases at low densities as found in white dwarf 
atmospheres, where conventional techniques fail to repro- 
duce gas opacities affected by non-ideal effects (e.g., the 
Lyman continuum opacity [30 . l5l| ) . The novel feature of 
the TSP formalism is that it allows one to evaluate, in 
a thermodynamically self-consistent way, populations of 
atoms and molecules in internal energy states as a func- 
tion of the available volume v, a variable linked directly 
with different degrees of interparticle perturbation. 

For one-component fluids, the TSP formalism is de- 
fined in terms of a Helmholtz free-energy model A based 
on the occupation number distribution N v of the variable 
v 133], 



A = 



N v . 
—r In 



N v \ l 
N 



dv + 



(3.6) 
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where j3 = l/(kT), A = V ' h2 P I '(27rm) , D is the space 
dimension, k the Boltzmann constant, h the Planck con- 
stant, m the particle mass and T the temperature. The 
first term on the right hand side of Eq. (|3.6j) is the ideal 
free energy as given by the TSP formalism. The space 
distribution N v for a particle system at equilibrium is de- 
termined by minimization of A subjected to total particle 
number and total volume constrains (see [33| for details) , 



N 



V 



N v dv, 



vN v dv. 



(3.7) 



(3.8) 



In Ref. [33j, the interparticle interactions in the fluid 
were incorporated into the TSP approach via the po- 
tential u v (w) and the PDF g v (uj). Therefore, the total 
potential energy $ was written in [33| as (|3.3|) . but as- 
suming that u v {uj) is independent of N v . In such case, 
$ is a linear functional of N v and the minimization free 
energy procedure yields 



N 

A 5 " 6XP 



1V + T + N ~ ^ 



with 



<t>v = n I u v (Lu)g v (uj)duj 



(3.9) 



(3.10) 



where the chemical potential, and 7 = P — $/V, the 
pressure minus the potential energy density, are the La- 
grange multipliers for our constrains. 

However, if we adopt the more general potential energy 
expression given by (|3.1[) . $ is a quadratic functional of 
N v . In this case, it is easily deduced that the minimiza- 
tion technique leads to 



N 
A 



D exp[-/3(7w + 0!, - fJ.)] , 



(3.11) 



with <p v and u v {u>) respectively given by (|3.10[> and p.2[) 
or, equivalcntly, 



n v 'U vv i (uj)g vv i (uj)dv'duj. 



(3.12) 



The main difference between Eqs. (|3.9[) and p. lip is the 
dependence of N v on <p v . It is important to note that the 
results obtained in [33| from the TSP approach for hard 
systems and van der Waals fluids remain unchanged if 
one uses (|3. 1 1 1) instead p.9[) . This is because of the spe- 
cial form of the two-particle potentials for these systems. 
However, Eq. p.lip is the correct one and must be used 
in general. In particular, Eq. p. lip is supported by 
results corresponding to attractive hard spheres (AHS). 
The implcntation of the TSP approach to the AHS model 



-i — 1 — 1 — r 



1 — 1 — 1 — r - 



H 1 1 I - 



r=i, x=i.05 
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1 1 1 1 
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0.1 



0.2 



0.3 



FIG. 2: (Color online) Contact values g(a + ) of the PDF for 
attractive hard spheres in the TSP formalism based on Eq. 
(|3.1ip (solid lines), are compared with exact values in the 
limit of zero density (triangles) , Monte Carlo simulations [H^l 
(circles), results of a perturbative theory (dashed lines), 
and results derived from Eq. (|3.9p (dotted lines). Here T* = 
kT/t is the reduced temperature, p = nd z the reduced number 
density, and A the potential range in units of the particle 
diameter d [i.e., £ = a(A 3 - 1) and a = 47rd 3 /3 in Eq. fDl))]. 



is shown in Appendix [D] We obtain the following contact 
value of the conventional PDF for this model, 



1 +\ Pi 
g{u} = a + = — < 

nB 



with /?7 and B given by Eqs. (|D5|1 and (ID6 
expansion of p,13p yields 



(3.13) 
A density 



g(a+) = e^ {l + [b + (2 + pe - 2e^)£\ n + 0(n 2 )} , 

(3.14) 

in agreement with the value expected at zero density 
g(a + ) = e -^"( a+ ) = e l3e . On the contrary, it may be 
deduced that the contact value at zero density result- 
ing from Eq. p.9p is incorrectly e^ e / 2 . Although the 
TSP approach is being developed mainly for the low den- 
sities found in stellar atmospheres, Fig. [5] shows that 
Eq. p,13p can also reproduce satisfactorily the results at 
moderately low densities obtained from Monte Carlo sim- 
ulations [52| . In conclusion, the correct expression p. lip 
for the space distribution N v at equilibrium in the TSP 
formalism, is reached with the minimization of a free en- 
ergy model based on the quadratic functional p. II) for 
the system potential energy and the use of the CPDF 
flw(w). 



IV. POISSON ENSEMBLES 

In this section we deduce an exact expression of g vv i (cj) 
for an infinite system composed of independent and ran- 
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t t I 

£dv vQv' £dv" 



FIG. 3: (Color online) Pairs (v,v') showing various geometri- 
cal quantities involved in the evaluation of the CPDF g vv i (u>) 
for Poisson ensembles in one and two dimensions. 



call Sp the probability for the occurrence of a pair (v,v') 
with the second point (say P'), which has an AV between 
v' and v' + Sv' , lying between the surfaces of volumes ui 
and uj + Slu centered in the reference point (say P), which 
has an AV between v and v + Sv. By definition of g vv > (w), 
Eq. (|2.5p . it is evident that the probability Sp can be 
expressed as follows 

Sp = n~ 1 n v n v >g l v ' J v }(uj)SvSv'5ui. (4.4) 

On the other hand, in the case of configurations class 
(iv) and according to (|4.2I) and (|4.3|) . the probability 
Sp equals the product of the following four probabilities 
(see Fig. [3]): 

(a) the probability (nSui) that the point P' (P) 
lies between the surfaces of to and uj + Suj centered in P 

(b) the probability (n^ v Sv) of finding a point be- 
tween the surfaces of v and v + Sv centered at P, where 

is the fraction of Sv that is not shared with v' , 



domly distributed points, the so-called Poisson ensemble. 
This system is equivalent to the thermodynamic limit of 
a fluid with non-interacting particles (the perfect or ideal 
gas, sometimes called Poisson fluid in the mathematics 
literature). Poisson ensemble serves as a prototype of 
more general point systems. 

The nearest neighbor distribution for random configu- 
rations of non-interacting points in three dimension was 
derived by Hertz [23[ . In the present notation, the distri- 
bution n v for Poisson ensembles in arbitrary dimension 
is given by 



n 2 e 



(4.1) 



This distribution can also be derived from the TSP 
method Indeed, (|4.1[) is the minimum Hclmholtz 

energy solution [Eq. (|3.11[) with cf> v = 0] for the distri- 
bution of available volume with the constrains of particle 
conservation (|3.7p and space normalization (|3.8p in the 
thermodynamic limit. 

For the evaluation of g vv '(uS), we use two basic prop- 
erties which characterize to a Poisson point process in 
arbitrary dimension and with average density n: 

(A) the probability that a given volume Q is found 

empty of points, except by a prescribed and 

finite number of them, is e~ nfl , (4.2) 



(B) the probability of finding a point within an 
infinitesimal volume dfl is ndfl. (4-3) 

We now apply these properties in the evaluation of the 
occurrence probability of pairs (v,v') in classes (i-iv). We 



(c) the probability (n£ v 'Sv') of finding a point between 
the surfaces of v' and v' + Sv' centered at P', £„/ being 
the fraction of Sv' that is not shared with v, 



1 - 



Sv' n v 



Sv' 



(4.6) 



(d) the probability (e nvUv ) that no points, apart from 
P and P' , lie in the union volume v U v' . 

For notational simplicity, explicit u> and v' (v) de- 
pendence has been omitted for (£„')• I n (&) an d ( c ) we 
have used the fact that, due to the meaning of available 
volume, there is no point into portions of Sv and Sv' 
which overlap with v' and v, respectively. Thus, for 
pairs in class (iv) we obtain 



Sp 



(n^ v Sv)(nCSv')(nSuj)e{uj - «+), (4.7) 



where Q(x) is the usual Hcaviside step function, which 
here takes into account that pair configurations with uj < 
max{j),w'} are forbidden [Eq. (|2.15[) ]. Therefore, using 
P~Tj) and that v Uv' = v + v' — v R v', it follows from (|4.4p 
and (JH7D that 



(™)/ \ 



0(w 



(4.8) 



A similar procedure applies inmediately to pairs (v,v') 
in classes (i-iii). In the case (iii) we must replace the 
assumptions (b) and (c) by the probability (nSv n Sv') 
of finding the shared nearest neighbor inside of intersec- 
tion volume (Sv PI Sv') of the spherical shells Sv and Sv' 
surrounding P and P' , respectively (Fig. |4j. Then 
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(4.9) 



7 



Area=_£a_6b_ 
sin(0+0) 




+ (^ + Uv')e(cj-v + )]. (4.14) 

The cross-section v <~) v' between v and v' for one- 
dimensional pairs (v, v') is easily deduced to be (v + v' — 
u>)/2 at w < v + v' and zero otherwise. Results for D = 2 
and D = 3 can be also deduced by simple geometry. 
Specifically, one finds that [HH 



/ l 



(v + v' — w), 
£[0-i S in(20)] 
"H!)'^ +^ [0 - | sin(20)l , 
|(2- 3 cos + cos s 9) 
k + ^(2-3cos<?(>-|-cos 3 



(D = l) 

(£> = 2) (4.15) 
(D = 3) 



for u + < w < cj^u' and v PI u' = for w > , where the 
angles and <fi are given by (D > 2) 



cos 



,2/D , ,,2/D 



2(w«)Vf 



(4.16) 



FIG. 4: (Color online) Shared neighbor configuration for D > 
2. 



where 7/ wu ' represents the ratio between the volume of 
intersection and the volume product of both shells, 



Sv n Sv' 

SvSv' 



(4.10) 



(the ui dependence of rj vv > is omitted for brevity). For 
pairs of class (i), i.e. mutual neighbors, the probabilities 
(b) and (c) are just considered by (a) and must be re- 
placed by Dirac delta distributions as expressed by Eq. 
(J33DJ). Th en 



9^) = \s{v-v')S{u-v)e^\ 



(4.11) 



where £ is the fraction of the overlap volume (v D v' /v) 
of two spheres with equivalent size and whose centers are 
separated by a distance equal to their radii (notice that 
£ depends exclusively on the dimensionality D). Some £ 
values are given in Table Hi] below. Similarly, in the case 
of pairs of class (ii), the probability of finding a nearest 
neighbor as given by either (b) or (c) is just included in 
the probability (a) and, according to (|2. 12|) . we have 



Hi) l \ ?t \ nvtlv' 



with 



V- 



min(u, v ). 



(4.12) 



(4.13) 



j2/D + V I2/D _ V 2/D 



ni/D 



(4.17) 



The quantities £„, £ v > r) vv i in ()4.14p can be evaluated in 
terms of elementary functions. For uj > u> vv ', with uj vv > 
given by Eq. (|2.13[) . we have = fu' = 1, r) vv ' = 0, 
v n v' = 0, and therefore 



9w'{u) = 9vv>(u) = 1, w > w vv >. (4.18) 
it can be shown easily (with help of 



For v + < to < ui vv i, 





and 



(D = l) 
(D = 2) 
(D = 3) 
(4.19) 



\5{uj — v — v'), 
(2 7 rv / Wsin[0 + 0])^ 1 , 
sinA^/V 2 / 3 sin(<9 + 



(D = l) 
(D = 2) 
(D = 3) 

(4.20) 

Some guidelines for the derivation of r] vv i are given in 
Appendix [A] Although r\ vv i was written in a particular 
way for D = 3, it is actually symmetric with respect to v 
and v' . Notice also that the prefactor 1/2 in (|4.19p and 
(|4.20p for D = 1 is related to the "edge surface" of the 
one-dimensional "volumes" v and v' (see Fig. [3]). 



Substitution of expressions (|4.8j) . (|4.9p . (|4.12l) and 
(|4.1ip into (|2.9p yields the CPDF for Poisson ensembles 



sw (w) 



S(v-v') £v_ . 
5 + ) S ( W - v +) 



A. Explicit expressions for D = 1 and D — > oo 



Substantial simplifications in Eq. (|4.14p occur in the 
limits of one and infinitely many space dimensions. In 
particular, the combination of (|4. 14[) with results (|4.15l) . 
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FIG. 5: (Color online) Evaluations of g vv i{uj) at dimensions 
D = l,2,3,oo for (nv,nv') = (0.2,0.2) (dotted lines) and 
(nv,nv') — (1,0.5) (solid lines). 
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FIG. 6: (Color online) Contributions g^f } (cj) and g (iv) {uj) 
(dashed and solid lines, respectively) to g vv i(uj) at v+ < uj < 
uj vv i , for D — 1,2,3, oo, and for the two (v,v r ) conditions 
shown in Fig. [S] 



(|4~T91 and (|T2TJ|) at D = 1, gives the CPDF for one 
dimensional Poisson systems 



e (nv)/2 e (™-)/2 

5— 8{v - v') + 



2n 



n{v+v' — u)/2 



5{uj — 7J + ) 



-Q(v + v' - lo)Q(lo - v + ) 

+ —5(lj-v-v') + Q(uj-v-v'). (4.21) 
2n 



On the other hand, as the spatial dimension tends to 
infinite the intersection volume of two off-center hyper- 
spheres («H u'), the overlap space between two off-center 
hypcrspherical shells (SvD 5v'), and the cross section be- 
tween a hypcrspherical shell and a hypcrsphere (Sv n v') 
no centered, all go to zero. Consequently, = — > 1 
and r] VV ' — > everywhere as I? — ?> oo, and therefore (|4.14p 
reduces to 



g vv >(u) = 



v') + - 

n 



5(u) — u+) + Q(oj — u+). 

(4.22) 

One of the main differences between the CPDF of Poisson 
ensembles at one and infinite space dimensions concerns 
the shared neighbor contribution. For D = 1, g™,(uj) is 
a singular function that corresponds to a single point of 
contact between v and v' , while that at D = oo the prob- 
ability density of finding a pair sharing the first neighbor 
tends to zero everywhere because r\ vv i — > 0. 



B. Discussion 

In Fig. Owe depict evaluations of g vv i (uj) based on Eq. 
(|4"7li|) for D = 1,2, 3, and oo, at (nv, nv') = (0.2, 0.2) and 
(nv,nv') = (1,0.5) (dotted and solid lines, respectively). 
For the same conditions, Fig. [S] shows the contributions 
to g vv '(u)) origined by pairs (v,v') in classes (iii) and 
(iv), as given by Eqs. (|4.9[) and (|4.8[) . The contribution 
which shows the largest variations is g^\uj), a direct 
consequence of the complex dependence of the overlap 
shell volume r\ vv i with the parameters v, v' and uj. 

As observed in Fig. [SJ g vv >(oj) vanishes for all di- 
mension at w < v+, which means that there are not 
very close pairs (v,v') as a consequence of the exclusion 
zones imposed by the AV v and v'. The main strik- 
ing features on the plots are the strong correlations be- 
tween pairs (v,v') at oj = u+ for all dimensions and at 
uj = uj vv i for D = 1 and 2. For, respectively, D = 1,2 
and 3, u> vv > = 0.4, 0.8 and 1.6 in the case of pairs 
(to;, Tit;') = (0.2,0.2), and aw = 1-5, 2.914 and 5.771 
in the case (nv, nv') = (1, 0.5). 

The singularity at uj = v+ , denoted by vertical lines in 
Fig. O is produced by pairs where at least one of the two 
points is the nearest neighbor of the other one, i.e., it has 
origin in pairs of type (ii) , together with pairs of type (i) 
for the example with v — v' (mutual nearest neighbors). 
The second strong correlation originated at one and two 
dimensions is due to the shared nearest neighbor configu- 
ration, i.e., pairs in class (iii). For D = 1, the singularity 
introduced by the 6(uj—v—v') function in Eq. (|4.2ip . seen 
as vertical lines in Fig. [51 is a consequence of the fact that 
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the shared neighbor must be located in the intersection of 
the outer edges of v and v', which for D = 1 corresponds 
to a single point in space. In the case of two-dimensional 
Poisson ensembles, the scaled overlap volume r] vv > gives 
nonzero values in the range v+ < uj < u) vv i and has a 
singularity at us = uj vv i , which is reflected in a dramatic 
rise of g m <(uj) (Fig. D = 2). 

At D = 3, the overlaping 77^^/ of the outer edges of 
i> and 1/ is less sensitive to variation of ui and therefore 
the contribution gjffi (uj) remains bounded for all u> (Fig. 

In this dimension, the discontinuity in g vv '(ui) at 
uj = uj vv i (Fig. [5]) is yielded by the sudden fall of the 
occurrence of shared neighbor pairs. For D = 00, the 
contribution of shared nearest neighbor configurations is 
everywhere negligible (dashed line in Fig. [5]) and that 
from pairs of class (iv) is uniform and equals unity (solid 
line), the total CPDF resulting flat at uj > v + (Fig. |5|). 



V. THE CPDF g v (uj) AT POISSON ENSEMBLES 

An expression for g v (oj) corresponding to Poisson en- 
sembles can be obtained in a similar, but simpler and 
more straightforward way than that previously used for 
g VV '(uj) (its derivation is left to the reader). This CPDF 
was previously introduced and analyzed in [33[ and reads 



g v (u>) = — S(uj — v) + 0(lu — v). 
n 



(5.1) 



Here, the Dirac delta function accounts for the nearest 
neighbor contribution of a reference point with AV v, and 
the Heaviside function describes the probability density 
of finding points (other than the nearest one) lying in 
the surface of uj > v centered in the reference one. It 
should be noted that (|5.1[) could be derived from (|2.7[) 
and (|4.14[) , but this procedure involves some analytically 
intractable integrals for D > 2. 

On the other hand, substitution of (|5.1|) into (|2.6I) pro- 
vides the expected PDF for Poisson ensembles, i.e., 



1. 



(5.2) 



This result reflects the well-known fact that no correla- 
tion exists between pairs of arbitrary points at Poisson 
ensembles. Clearly the departures from unity in g v (uj) 
and g vv > (ui) are due to that the available volumes v and 
v' introduce correlations between pairs at short distances, 
which correspond to u < v and uj < uj vv i for g v (uj) and 
flw(w), respectively. 



A. Relative importance of the different 
contributions to the CPDF g v (u>) 

The results obtained in the derivation of g vv ' (uj) can be 
used to evaluate various statistical neighboring functions. 
In particular, expressions (|2.6p and (|2.7[) can be applied 



to each class (v = i, ii, hi, iv) of pairs (v, v'), 



g<r>(cj) = - n v9 M(uj)dv. 



(5.3) 



(5.4) 



The quantities in Eqs. (|5 .3[) and (|5.4[) express partial 
probability densities that describe the spatial distribu- 
tions of specific pair configurations. For mutual neighbor 
pairs, from (|4.11[) . ()5.3|) . and ()5.4[) we obtain 



g®(w) = ±e«- 1 ~> n "6(u>-v), 



g^(uj) = e^ nu . 



(5.5) 



(5.6) 



Here, it is useful to separate pairs in class (ii) into sub- 
classes (iia) and (iib). Pairs (v,v') where the point with 
AV v' is the (non mutual) nearest neighbor of the point 
with AV v belong to the subclass (iia), otherwise they 
belong to the subclass (iio). Then, Eq. (|5.3p yields 



gM(w)=g!r>{u>)+g}r>(u>) 



(iia) 



(iib) ( 



(5.7) 



with 



and 



ti™\uj) = 5{uj - v) / &e- nC »'dv', (5. 
Jo 



gi iib \u)=[£ v e- nC «'].e(u;-v), 



(5.9) 



where C v ' = v' — v PI v' is the volume of a crescent in 
the intersection of v and v 1 . The term g^ la \ui) is the 
probability density of finding a non-mutual, first neighbor 
in the surface of ui centered in a reference point with AV 
v. The term gi (u>) is the contribution to g v (uj) due 
to pairs where the reference point with AV v is the non- 
mutual nearest neighbor of the other one. Note that the 
Dirac delta term in (|5.ip arises from contributions due 

to g£\ui) and gi" a \uj) through (|2.7I) . which then implies 
that the integral in (|5.8j) is given by 



Z v >e~ nC *'dv' = - 



1 _ e (C-l)r 



(5.10) 



The contribution to the PDF g(ui) due to pairs in sub- 
classes (iia) and (iib) may be found from (|5.4I) . (|5.8p and 
(|5.9p . Thus, using (|5.10p . one obtains 



g [iia) {uj) = e~ nw - e ^- 2)nu . 



(5.11) 



By symmetry, pairs in subclasses (iia) and (iio) have 
equivalent PDFs, i.e., g( nb \ui) = g(" a )(uj), which implies 



(uj) = I 




-nvUv 



dv = e~ nu _ e (c - 2)n ". 



(5.12) 
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FIG. 7: (Color online) Contributions due to (ii&), (iii) and 
(iv) pair configurations (solid, short dashed and long dashed 
lines, respectively) to the CPDF f/„(w) as a function of v, for 
one dimensional Poisson ensembles and different distances as 
measured by w (values raw indicated on the plots). The total 
CPDF is shown by the dotted line. 
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FIG. 8: (Color online) Same as Fig. [7] for D = 2. 



A similar analysis applies to pairs in classes (iii) and 
(iv) but their contributions to g v (u>) and g(oj) cannot be 
determined analytically except for D = 1 and D = oo. In 
general, these functions are expressed as integral forms, 

gM(w) = r] vv ,e- nC ^dv'e(cj - v), (5.13) 
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g« v) (w)=n Uve-^'dv'Qiu-v), (5.14) 
Jo 



g^ u \uj)=n / Vvv ,e- nvUv 'dv'dv, (5.15) 
Jo Jo 



&tv>e- nvUv dv'dv. (5.16) 



Again, by virtue of (|2.7j) . the Heaviside term in (|5.1j) 
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comes from ([SU]) . (f5TT3|) and (f5~T4"l) . therefore 

fl c«*)( w ) + = [i - (e,e-" c -')^=J e(« ~ «)■ 

(5.17) 

which implies 

/" (VW + nUv>)e- nC *'dv' = 1 - [f,e-" c .'] , ^ . 
Jo 

(5.18) 

Similarly, with (f5TT2|) and (|5~T7| we obtain 

^(w) + = 1 - 2 e -" w + e^" 2 ^". (5.19) 

Explicit mathematical formulas for gi (w) and (uj) 
functions at D = 1 and D = oo have been collected in 
Appendixes [B] and [Cj 

The explicit expression for g"> (uj) at D = 3 [as given 
by Eq. (|5.6[) ] was previously obtained by Larsen & Stratt 
|29|. To the best of our knowledge, all other results de- 
rived here are new. 



B. Discussion 

Figures [5J [HI and [TO] depict the partial probability 
densities g^ lh \uj) (solid lines), gi m \ui) (dashed lines) 

and gi m \uj) (long dashed lines) as functions of the AV 
v (0 < v < uj), for different volume uj and dimensions 
D. These evaluations come from Eqs. (|5.9p . (|5.13|) . and 
(|5.14[) . It is worth to notice that the sum of these three 
contributions gives the unity (dotted lines in the figures), 
which is actually the value of the total CPDF g v (oo) for 
v <lo [Eq. ([5~T]) ]. 

At finite dimension and short pair separations (nuj <C 
1), the CPDF g v (u>) is almost entirely composed of con- 
tributions due to pairs in classes (ii&) and (hi), i.e., the 
reference point with AV v is the nearest neighbor of the 
other point at uj or both have a shared nearest neighbor. 
In the one dimensional system (Fig. UJ), shared neighbor 
configurations (short dashed lines) are somewhat more 
likely than those in class (ii&), while for D > 1 the ma- 
jor contribution arises from pair of class (ii&) (solid lines, 
Figs. [3] and [9]) • Pahs sharing the nearest neighbor be- 
comes highly unlikely at high space dimensionalities for 
any pair separation (Fig. [TO)) . As expected, class (iv) 
(long dashed lines) is the dominant pair configuration at 
any D for great distances (nuj > 1) from the reference 
point. 

At finite dimension, as the available volume v of the 
reference point increases with fixed uj, the surface frac- 
tion in Eq. (|5.9p increases and the exclusion volume 
C v i =tAl decreases. Consequently, the probability of finding- 
pairs at class (ii&) increases monotonically with v. On the 
contrary, the chances of finding four-point configurations 
(pairs at class iv) reduce with increasing v (or decreasing 
w) because the surface fraction at the integrand of 
Eq. (|5 . 14[) takes on average lower values (which are not 
fully compensated by the increasing of the other terms in 




FIG. 11: (Color online) Contributions g^ v \oj) due to v =i, ii, 
iii and iv pair configurations (solid, short dashed, long dashed 
and point-dashed lines, respectively) to the PDF g(uj) vs the 
volume uj, for Poisson ensembles at dimensions D = 1,2,3 
and oo. The total PDF is shown by a dotted line. 



the integrand). The higher the Euclidean space dimen- 
sionality, the weaker the v dependence of the geometrical 
quantities (£„, rj vv t , C v >), and the partial probability den- 
sities gi^ (uj) analyzed here become independent of v at 
the limit D oo (Fig. [TO]) . 



C. Relative importance of the different 
contributions to the PDF g(uj) = 1 

Figure [TT1 illustrates the behavior of the partial proba- 
bility densities g^\uj) given by Eqs. ([H]), (|5TTT]) . (pTTOj) . 
(j5"T5l) and (j5~TO|) . Note that g {ii \ui) is given by the sum 
of (|5.1ip and (|5.12[) . The term g^{u) dominates the 
regime of low volume (nuj <C 1) at any D, i.e., most 
points with a first neighbor at very short distance (lower 
than the mean separation between points in the ensem- 
ble) form pairs of mutual nearest neighbors. At interme- 
diate distances from a reference point (which correspond 
to nuj « 1), there is a relatively narrow uj range where 
pairs in class (ii) are the main contribution to g(w) for 
any D. As the distance from the reference point increases 
(nuj > 1), pairs in class (iv) become the dominant con- 
tribution to g(uj). Indeed, g^ lv \uj) is a monotonically in- 
creasing function of uj and tends asymptotically to unity, 
while the remaining contributions go to zero for uj large 
enough. The term g( IM ) (oj) corresponding to pairs with 
the shared nearest neighbor is never dominant and has a 
maximum contribution of 18.39%, 16.72% and 14.54% 
at nuj = 1, 1.226 and 1.370 for D = 1,2 and 3, re- 
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TABLE II: Average fractions of pairs in various classes 



D 


a 


C 


f{i) 


j(iia) j(iib) 






1 


2 


1/2 


0.667 


0.333 


0.500 


0.50 


2 


7T 


4,r-3 3 / 2 

5/16 


0.622 


0.378 


0.618 


0.63 


3 


4tt/3 


0.593 


0.407 


0.709 


0.71 


CO 








0.500 


0.500 




1.00 



* Values listed in Ref. [58| 



TABLE 111: Mean cj-volumes and average distances between 
points of pairs in classes (i), (ii) and (iii). (r) and (uj) are 
expressed in units of mean distance between points n -1 ' 25 
and mean volume per point n , respectively. 



D 


{r) d) 


(r) {u) 


/ r \(«i) 








1 


0.333 


0.833 


1.000 


0.667 


1.667 


2.000 


2 


0.394 


0.674 


0.880 


0.622 


1.622 


2.797 


:', 


0.466 


0.683 


0.909 


0.593 


1.593 


3.961 



spectively. This maximum value gradually decreases and 
shifts to larger uj with increasing dimension D. The con- 
verse behavior is followed by g^(uj). Indeed, the max- 
imum probability for the occurrence of pairs where one 
and only one of the two points is the nearest neighbor of 
the other occurs at nuj = 0.811 (29.63%), 0.781 (34.67%), 
0.761 (38.06%), and 0.694 (50%) for D = 1, 2, 3 and oo. 



D. Pair fractions and mean separations 

Finally we mention that one can use g vv r(uj) to obtain 
other useful measures of the structure of Poisson ensem- 
bles or many-body systems in general. For examples, the 
fraction f' v ) of points in classes v =i, ii and iii, averaged 
over the whole ensemble, and the mean pair separation 
(r)( v > in each class. They are defined as 



f [v) = / ng^\u)du, 
Jo 



(5.20) 



and 



(r) lv) = j& J nr(L>)gM(u,)du, (5.21) 

with r(uj) = [uj/a) 1 / . Results for pairs in classes (i) 
and (ii) at Poisson ensembles are found analytically using 
Eqs. (EH), H5HU) and ([BTT2) . Thus 



/ 



2-C 



j(Ha) j'(iib) 



i-C 



«(2-C)] 



l/D ' 



(r 



w = Ell±i) [(2-0 1+1 / g -i] 

[na(2-C)] 1/D (1-C) 



(5.22) 



(5.23) 



(5.24) 



(5.25) 



Notice that the nearest neighbor of any point is cither 
mutual or non mutual one and therefore /« +/(" a ) = 1. 
Mean values {u}^"' can be also computed in similar form 
to (r>M, 



.(2-C)' 



, (ii) 



(5.26) 



Clearly, results based on the volume uj are more concise 
than those upon the distance r. 

The fraction / W of mutual nearest neighbors at Pois- 
son ensembles was initially studied by Clark [55|], al- 
though the exact evaluation [equivalent to that given by 
Eq. (|5.22p 1 was obtained by Dacey [56| and Cox [57|. The 
fraction of pairs sharing the nearest neighbor has 

also been computed by Schilling j58j. All these authors 
have used techniques different from that adopted in the 
present work, which are based on conditional pair dis- 
tribution functions. Our evaluations of f^ v \ (rp v > and 
{luY v \ together with values of the parameter a and the 
overlap volume £ are listed in Tables [TT| and IIIII Results 
corresponding to pairs in class (iii) were evaluated by nu- 
merical integration of (|5.20[) and (|5.21j) . The fractions 
f( m ) computed here are in agreement with those ob- 
tained by Schilling [58[ (see Table ITT|) . Interesting enough, 
Schilling showed that tends to unity as D — >• oo. 

This implies that whereas g^™) (to) tends to zero every- 
where (see Section TlV Ap with D — > oo, its integral over 
uj [Eq. (|5.20p ] converges to unity. 



VI. 



CONCLUSIONS 



A new mathematical descriptor of spatial structures in 
many-body systems have been proposed. Specifically, we 
have introduced a two-point distribution function, de- 
noted by g vv i(u>) and defined by Eq. (|2.5p . which de- 
scribes the probability density of finding pairs of objects 
with each object having a nearest neighbor at a cer- 
tain distance. For simplicity in the mathematical frame- 
work, distances between objects have been substituted 
by spherical volumes whose radii denote the separation 
between two objets of a pair (volume w), and between 
each object and its nearest neighbor (volumes v and v'). 

The conditional pair distribution function (CPDF) 
g vv r (uj) reveals a richer structural information than that 
provided by the conventional PDF g(u>). Indeed, the vol- 
umes v and v' introduce correlations in the CPDF which 
are not implicit in g(uj). Moreover, nearest neighbor in- 
formation from v and v' let us to identify four types 
of pairs (mutual nearest neighbors and neighbor shar- 
ing pairs, among others), which provide a considerable 
amount of information concerning the local microscopic 
density and spatial relationships in many-body systems. 

General relations were established among g vv '(uj), the 
standard PDF g(uj), and g v (uj). The function g v (u), pre- 
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viously introduced in [33J , is a CPDF in which the near- 
est neighbor information is only considered for one ob- 
ject of each pair. We have shown that, using appropriate 
many-body interaction potentials, g v {^>) and g vv >(uj) can 
be used to evaluate the total potential energy of a fluid, 
which reduces to that typically adopted in liquid theory 
with an averaged pairwise potential. It was shown also 
that the potential energy expression based on g vv >{uS) is 
the correct one to be used in the so-called thermostatis- 
tical space partition (TSP) formalism (33l. |3~H- 

Wc have studied the function g vv i(ui) for Poisson en- 
sembles and derived its exact expression in any space 
dimension. We have obtained explicit results for g vv > (uj) 
at D = 1, 2, 3, oo, and found simple close forms for both 
D = 1 and D = oo. To our knowledge, all these re- 
sults are new. Furthermore, functions g vv i(uj), g v (u>) 
and g(uj) corresponding to Poisson systems have been de- 
convolutcd into contributions from the four pair classes. 
With the exception of the function g^(ui) given in [29[ 
for mutual nearest neighbors at D = 3, all the partial 
probability densities obtained in the present study are 
new. Our theoretical results should have interesting and 
practical applications in the study of neigborhood struc- 
tures in a variety of areas, such as physical and biological 
sciences, and sociology. The extension of present evalua- 
tions to non-Poisson processes is open to future efforts. 



a distance r. This intersection volume is a torus with 
section area SaSb / sm(9 + <f>) and perimeter 2ttR, where 
9 and <j) are defined as before and R = a sin 9 = b sin </>. 
Therefore, we can write 



Sv n Sv' 



2na sin #(5a(5 



CD = 3). 



(A2) 



Equation (|420|) at D = 2 and D = 3 results from (f4~10|) . 
(|A1[) and (|A2|) . once the radii (a and b) and thickness 
(Sa and Sb) are written in terms of v = aaP , v' = ab D , 
Sv = Daa D ~ 1 8a and Sv' = Dab D ~ 1 Sb, with a given by 

(ED. 



Appendix B: Partial probability densities at D = 1 



For one-dimensional Poisson ensembles (D = 1), Eqs. 
(j4~2l"T) . and (JEI1) yield 



Acknowledgments 

The authors are indebted to Professor Andres San- 
tos for helpful suggestions and comments. This work 
has been supported by CONICET (Argentina) through 
Grant PIP 112-200801-01474. 



Appendix A: Covolume of two shells 

In this Appendix we give a brief derivation of the scaled 
overlap volume r) vv i of shells Sv and Sv' surrounding the 
available volumes v and v' [Eq. (|4.20[) ]. This quantity is 
directly related to the shared neighbor configuration in 
pairs (v, v'). The shared nearest neighbor of a pair (v, v') 
must be located in the intersection of the outer edges 
of v and v' . For D = 1, this intersection corresponds 
to a single point in space, which yields the Dirac delta 
function at (|4.20[) with a prcfactor 1/2 because the shared 
neighbor is located in one of the two sides of v. At D = 2, 
Sv H Sv' is the overlap surface of two off-center circular 
rings with thickness Sa and Sb, as shown in Fig. 0J It is 
shown easily by simple geometry that 



Sv n Sv' 



2SaSb 



sin 



(D = 2), 



(Al) 



where 9 and (j> are defined in Fig. U and expressed by 
Eqs (|4TTo| and (|4TT7|) . In the case D = 3, 5v n Sv' is 
the overlap volume of two spherical shells with internal 
radii a and b, thickness Sa and Sb, and centers separated 





1 

= — e 

n 


~ nv ' 2 5{oj - v), 




g ( : ia) {u) 


-k( 


1 - e -'"' /2 ) 5(lo - 


v), 


g ( : ib) n 


1 

2 6 


- nuj+nv / 2 e(uj - v), 




gf'Hu) 


1 

2 6 


-«(<--«)e( w - v), 






= 1 


- i (er^-^ + e 




and 









(Bl) 
(B2) 
(B3) 
(B4) 

(B5) 





H 


— 3nu/2 
e ) 


g(iia) 


H 


= 5 ( " fc) M 


g(iii) 


H 


- — c-™ 
" 2 




M 





- 2e- nu + e 



-3nu /2 



(B6) 
(B7) 
(B8) 
(B9) 



Clearly, the sum of terms g^ (u) and g^ (lu) over all pair 
classes (y =i, iia, ii&, iii and iv) let us recover the exact 
expressions for g v (oj) and g(oj), Eqs. (|5.ip and (|5.2[) . 
respectively. 
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Appendix C: Partial probability densities at D — oo 

For Poisson ensembles at infinity dimensions, Eqs. 
(l4~2H . and (USD yield 



and 



<? ( " a) M 
<? ( " l) M 



(1 - e"™) *(w - w), 



e-" w 6(w-i;), 



(l-e-™)e(u;-v), 



-2nuj 



(iib) 



( W ) 



1 — 2e _naJ + e~ 2naj . 



-2nu 



(Cl) 

(C2) 

(C3) 
(C4) 
(C5) 



(C6) 
(C7) 
(C8) 
(C9) 



As in the case D = 1 (Appendix |B|) . Eqs. (|5.1[) and ([5 
are recovered from (|CHIC5[) and (|C6tiC9|) , respectively. 



Appendix D: Attractive hard spheres in the TSP 
approach 

The model of attractive hard spheres can be used to 
check the role played by the CPDF in the evaluation 
of the potential energy contribution within the TSP ap- 
proach. We consider a pair interaction potential which 
consists of a hard core repulsion together with an attrac- 
tive square well, 



with a, £ and e constants. For the present purpose, we 
may use the expression of g v {w) given in |33j [eq. (41)] 
at the low density limit, i.e., 

g v (u) = -5(w-v-a*) + Q(u)-v-a*), (D2) 
n 

where 0(x) is the Heavisidc step function and a* = a/2 
is a reduction of the AV of hard particles in the limit of 
low density. Thus, with ([Dl) into Eq. (|3~T0|) we obtain 

-co, v < b, 

-e[l + n(b + £ - v)] , b<v<b + £, (D3) 
0, uj>b + t 



where b = a - a* . With help of Eqs. and (|3T8|) 

in the thermodynamic limit (N, V — > oo with N/V = n 
constant) one finds 



and 



with 



n/3j 

( 

B 



P[j(v-b)+<j>„] 



Pi = n-£, 



D 



/3e(l+n£) 



1 _ e -/8(7+«e)£ 
Pil + ne) 



(D4) 
(D5) 

(D6) 



and 



C = [l + p 7 (b + 0]e™ + 

x {1 + j8( 7 + ne)b - [1 + 0(7 + ne)(6 + f)] 



/3 7 



/3e(l+nJ) 



x e 



-/3( 7 +ne)£ 



} 



(D7) 



oo, uj < a, 
u(uj) = ^ — e, a < uj < a + £, 
0, w>a + ^, 



The contact value g(w = a + ) of the PDF given in Eq. 
(Dl) ()3.13[) is directly obtained from Eq. 
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